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We study the physical content of the electromagnetic vacuum energy of a random medium made of 
atomic electric dipoles. First, we evaluate the contribution of statistical fluctuations to the average 
total vacuum energy, which is made out of the integration of the variations of the Lamb shift with 
respect to the coupling constant. While the Lamb shift is a function of the electrical susceptibility 
only, the vacuum energy is generally not. Second, we make clear why the effective medium bulk 
energy does not account for the total vacuum energy of a molecular dielectric. Consequently, the 
Lamb shift does not derive from the effective medium bulk energy except at leading order in the 
molecular density. The local field factors provide natural cutoffs for the spectrum of the total 
vacuum energy at a wavelength of the order of the correlation length. Third, we investigate to what 
extent shifts in the spectrum of the dielectric constant may be attributed to the binding energy 
of a dielectric. In particular, in the continuum approximation we have found a relation between 
the electrostatic binding energy and the Lorentz-Lorenz shift. Nonetheless, we conclude that the 
knowledge of the spectrum of the refractive index is insufficient either to quantify the energy of 
radiative modes or to estimate the electrostatic binding energy of molecular clusters. 



I. INTRODUCTION 

It has long since been recognized that the Casimir ef- 
fects, the van der Waals (vdW) forces and the Lamb 
shift share a common origin [l|-|3|. It is customary to 
ascribe the vdW forces and the Lamb shift to the short- 
range interactions between the microscopic constituents 
of dielectric media. In contrast, the Casimir effects are 
attributed to the long-range interactions between macro- 
scopic objects. Also, some authors refer to the energy of 
discrete modes as Casimir energy while they term bulk 
energy that of the continuous spectrum In this 

paper we deal with a translation-invariant molecular di- 
electric made of atomic electric dipoles. Therefore, aside 
possible geometrical resonances in clusters, the spectrum 
of electromagnetic (EM) fluctuations is continuous and 
the above distinctions are unnecessary. 

We aim to compute the EM energy of the ground state 
of the dielectric. It has been also recognized that this en- 
ergy can be interpreted in two alternative ways. That 
is, when speaking of the interaction energy of fluctuating 
dipoles it is referred to as binding energy 04^] . Alterna- 
tively, when computed out of the variations of the zero- 
point vacuum fluctuations of the EM field it is termed 
vacuum energy [Hi. Evidently, this distinction is 
just semantic. 

The EM energy manifests in the shifts of bound state 
levels and in the shifts of the spectrum of the dielectric 
constant. However, there does not exist an obvious cor- 
respondence between them. The binding energies of the 
mcchum prior to the coupling of the dipoles to radiation 
are those of the quantum atomic states and those of the 
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short range forces which structure the spatial configura- 
tion of the dipoles. The former can be parameterized by 
the bare resonant frequencies of the atomic transitions. 
The latter depend on the precise form of the interaction 
potential between the molecules and reflect on the spa- 
tial correlation functions which give rise to the spatial 
dispersion of the dielectric function. We will restrict our 
analysis to situations in which the correlation functions 
do not depend on dipole-dipole interactions. Therefore, 
the energy shifts will affect only the atomic levels, even 
though they will depend on external degrees of freedom 
(d.o.f.). For this reason we will refer to the energy shifts 
as generalized Lamb shifts. They are made of two con- 
tribution. The first one is common to all the dipoles and 
is due to the coupling of each isolated dipole to bare ra- 
diation. That is the free-space Lamb shift [ll| - [l3i , which 
is additive and whose integration is a self-energy termed 
free space Lamb energy. The second contribution is due 
to the interaction of each dipole with the rest, which is 
mediated by the multiple scattering of virtual photons. 
That is the scattering Lamb shift [12, H, 113 ■ Physically, 
it is the energy needed (or released) in the removal of one 
of the dipoles out of the dielectric. The vdW forces and 
torques acting upon each dipole derive from it [HI, [l^ . 
Its integration is a non-additive binding energy which in- 
volves a number of dipoles and collective d.o.f. The non- 
additive character of the scattering Lamb shift can be 
better explained in terms of the decay of excited atoms. 
When one of the dipoles is in an excited state, its fre- 
quency is renormalized by a 'resonant' Lamb shift which 
contains a scattering contribution [l^, [l^- There- 
fore, part of the energy released in its decay belongs to 
collective d.o.f. For this reason, when several neighbor- 
ing dipoles are excited and decay at a time the energy 
released is not the sum of their individual energy shifts. 
On the other hand, according to Schwinger's source 
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theory, the vacuum energy can be computed out of the 
integration of the variations of an effective action for the 
EM field w.r.t. some effective coupling to the dielectric 
medium [l7l - [l9| . Since the variations are taken adiabatic, 
the vacuum energy only depends on the initial and final 
properties of the dielectric. However, the result depends 
on both the action varied and the parameter of the vari- 
ations. 

In principle, there is no obvious correspondence be- 
tween the shifts on the dispersion relations of optical 
modes and the Lamb shifts. To this respect, Feynman 
[20j . Power [ll| and Milonni '21S have proved that the 
free-space Lamb shift can be computed from the vari- 
ation of the vacuum energy due to the presence of a 
single atomic dipole. That result has been extended by 
Schaden, Spruch and Zhou 22] to the case of a uniform 
distribution of non-niutually interacting dipoles. In refer- 
ence [3| Milonni, Schaden and Spruch have extrapolated 
those approaches to the computation of the Lamb shift 
from variations on the effective medium bulk energy of 
a molecular dielectric due to a change in the refractive 
index. However, that result was not totally conclusive 
since near field effects and local field factors (LFFs) were 
neglected in the calculation. In fact, the findings in 
suggest that that result does not hold when proper ac- 
count of LFFs and inherent correlations is taken. 

The paper is organized as follows. In Sec. |TT]we define 
fundamental concepts and review the basic formulas ob- 
tained in [23II. In Sec. Illll we explain the role of statistical 
fluctuations in the computation of the vacuum energy of a 
random medium. In Sec. llVl we explain why Schwinger's 
approach is not suitable to study the Lamb shift. In Sec. 
IVl we compute the Lamb shift and the vacuum energy 
density up to order two in the molecular density for a 
hard sphere model. The results are compared with the 
approach of Milonni et al. Q in Sec lVIl Sec. I VIII deals 
with the EM binding energy of a molecular dielectric in 
the continuum approximation. We analyze previous ap- 
proaches [1, H^l and comment on the possibility to find 
signatures of the vacuum energy on the spectrum of the 
dielectric constant. The conclusions are summarized in 

Sec. rvml 

Regarding notation, we will label three-spatial- 
component vectors with arrows and three-by-three ten- 
sors with overlines. We will denote the Fourier-transform 
of functionals with q, w-dependent arguments instead of 
the r, ^-dependent arguments of their space-time repre- 
sentation. Quantum operators will be denoted with bold 
letters and hats on top. 



II. GROUND STATES, POLARIZABILITIES 
AND PROPAGATORS 

A generic statistically homogeneous molecular dielec- 
tric is made of N equivalent point electric dipoles 
in a volume V, with a typical correlation length ^, 
such that, in the limit V ^ V ^ the aver- 



age numerical density is uniform, Jyd^rp{r) = 
p = N/V. For the sake of simplicity we will con- 
sider the canonical ensemble at temperature T with 
Za'(T) = /nLd'i?'exp[-C/(i?i,..,i?^)/fcBr] the N- 
body canonical partition function. U is the interaction 
energy, which is considered a function of the external 
d.o.f. of the dipoles only. That is, of the center of mass 
positions and velocities, Ri,Vi, i — 1,..^N. Eventually, 
in the limit ksT <^ miV^,U™"'^, rrii being the scatterer 
mass and the potential barrier from which steric 

forces derive, we can neglect the scatterer dynamics and 
assume that the potential U is short-ranged. This way, 
we can ignore Doppler effects and guarantee that the 
correlation functions do not depend on the long-range 
dipole-dipole interaction which is the object of our study 
[Jll^l- The dielectric constant and rest of effective opti- 
cal parameters are stochastic functions which admit clus- 
ter expansions. That is, series of n-scattering terms com- 
puted out of the convolution of n-body spatial correlation 
functions with single particle polarizabilities and electric 
field propagators -see below and [Hj. Further, we will 
work for simplicity in the zero-temperature limit, T = 0. 
In addition to the aforementioned inequalities it demands 
fcsT <C hujQ, with fvjjQ the typical excitation energy of 
the atomic dipoles, so that all the dipoles remain in their 
ground state. 

The interaction Hamiltonian in the electric dipole ap- 
proximation reduces to 

H^nt{t) - j ei{T,t) ■ i(f,i), (1) 

where e is the electronic charge an ef(r, i), E(r, t) are the 
dipole moment density and electric field operators respec- 
tively in the Heissemberg representation. Before statisti- 
cal averages being performed, we consider a fixed configu- 
ration of indistinguishable two-level atomic dipoles for 
which -ef(f,i) = -e^i{t)5'^^\r - R'), with e?,(t) 

the dipole moment operator associated to the dipole with 
position vector i?'. The EM radiation is intended as a 
reservoir of infinite number of d.o.f. while the internal 
d.o.f. of the dipoles constitute small systems. The cou- 
pling between the radiation reservoir and the dipoles is 
weak so that ordinary time-dependent perturbative cal- 
culations can be carried out. In the following we give 
some definitions and explain the notation. 

A. EM vacuum and atomic ground states 

Prior to turning on the interaction Hamiltonian, we 
denote the translation-invariant EM vacuum at T = 
by \Qo) and the bare ground and excited atomic states, 
common to all the dipoles, by \Al) and \Bl), i — 1, ..,N, 
respectively. Dipole and EM fluctuations are uncoupled 
and uncorrelated. 

Next, consider the dipoles infinitely separated from 
each other or otherwise a unique dipole in free space. 
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say the one at i?% and turn on the interaction Hamil- 

tonian —eri{t) ■ E(i?*,t). In this case the EM vacuum 
gets polarized locally only due to the dipole fluctuations. 
Reciprocally, the atomic states get renormalized by the 
EM fluctuations and radiation reaction, being denoted 
by \Aq), I-Bq), which are common to all the dipoles. 

Further, when the N dipoles are brought together to 
form a specific configuration, say m, with position vec- 
tors {Rln}, i = 1,..,N, the EM vacuum gets polarized 
by the local and non-local dipole fiuctuations and the 
atomic states get renormalized by both the local interac- 
tion of each dipole with the 'bare' EM fluctuations and 
the non-local mutual interactions between the fluctuat- 
ing dipoles. Alternatively, it can be interpreted that each 
atomic state gets renormalized by the local coupling of 
the corresponding dipole moment to the polarized fiuc- 
tuations of the EM field -i.e. polaritons. In this case the 
polarized vacuum is denoted by \ flm,), which is not trans- 
lation invariant, and the renormalized atomic states are 
denoted by \Al^), i — 1,..,N, being aU different 

in general. The states \flm), 1^™), l^m) '^iH be defined 
later on by their physical content in a way analogous to 
that in dl i^. That is, by the fluctuations of the EM 
field and dipole moment operators. 

Finally, consider a dielectric as a random medium de- 
scribed by a statistical ensemble of dipole configurations. 
As mentioned above, those configurations depend only 
on external d.o.f. and are independent of the atomic 
states. Now, on top of the quantum fluctuations that 
the dipoles and the EM field induce on each other, there 
are the stochastic fiuctuations induced by the random 
distribution of scatterers. However, generally stochastic 
fiuctuations do not enter the problem in the same fashion 
as quantum fluctuations. While the correlation time of 
the latter, r, satisfies r ^ ^^ab^ where hijj\g is the typ- 
ical energy transferred in the internal processes, so that 
quantum fluctuations drive the dynamics of the atoms; 
the correlation time of the ensemble of configurations is 
much longer. Thus, the situation is generally that for 
quenched disorder and stochastic fiuctuations enter only 
as performing ensemble averages -see also note 4. Let 
be the probability density of the mth dipole configuration 
at temperature T, 



N 



N 



with Z7v-i|^i^p(T) the conditional {N — l)-body parti- 
tion function, 



N 

•' i=i 



(4) 



The statistical average of the expectation value of any 
local operator O in the zero-temperature limit, = 
(Aj„,ri„|a|a„,AL)> reads. 



N 



P. 



■T=0\ 



(5) 



The choice of i is actually irrelevant since the dipoles 
are statistically equivalent and so must be any average 
observable. Alternatively, we can express symbolically 
i^)'avg' the expectation value of O in the stochastic 
polarized EM vacuum and renormalized atomic ground 
state, |riai,g, Aavg)- Formally, this can be written in terms 
of the stationary reduced local (i.e. at each dipole loca- 
tion) density matrix of the system. 



N 



P. 



T=0| 



J = l 



(6) 



in which the integration amounts to taking the partial 
trace over the statistical mixture of dipole configuration 
states. The averaged expectation value of O at T = is 



(O) 



avg \ J 



N 



J = l 

)|,-=Tr{p^=°(f).0}, 



Pi n d'^™ = '^n\T) n A^K. exp [-Uiiil, .., i?^)/fcBT] 

(2) 

and -Pml_Ri =F the conditional probability density with a 
dipole fixed at f -say the ith. 



N 
3 = 1 



N 



S^^Hr-Rl^)Pll[^'RL, (3) 



where \Al^,V,"^) above depends implicitly on {-R^}, j — 
I, ..,N. If O is a time-dependent interaction which cou- 
ples weakly the dipoles to the EM field, by writing p as a 
direct product of EM and atomic states we assume that 
the EM field behaves as a stationary reservoir w.r.t. the 
weak interaction and the time-correlation between the 
reservoir and the dipoles is negligible as considering the 
dynamics of each dipole -cf. Ch.IV of [131 ■ Also, in a 
statistically homogeneous dielectric the spatial correla- 
'tions functions are translation invariant and so are the 
stochastic states and average expectation values. Hence 
it suffices to compute them at a single point. 



B. Response functions and fluctuations 

Our goal is in the first place to give an expression 
for the Lamb shift of the atomic ground states, S'"^^. 
Second, we will use that expression to compute the total 
binding energy of the ground state of the dielectric at 
zero temperature. Because the latter can be equally 
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intended as a variation of the zero-point energy of the 
EM vacuum due to the presence of the dielectric -in 
the same spirit as the interpretation of Power [ill ]. 
Feynman ]20| and Milonni [2l[ for the case of an only 
atom in free space- we will refer to it as total vacuum 
energy density, . We will reinforce this interpretation 
giving an expression for in terms of the source field 
propagator, in accordance to Schwinger's formalism. 

In order to compute the Lamb shifts and the total 
vacuum energy it is necessary to know the expressions 
for the two-time quadratic correlation functions of the 
EM field and dipole operators in their renormalized 
vacuum and ground states respectively -cf. Ch.IV of |23] 
and [2^. Generically, once the interaction Hamiltonian 
is turned on, ordinary second order time-dependent 
perturbation theory [29| yields an energy shift in the 
ground state of each dipole, say the ith one with position 
vector i?*, which corresponds to the generalized Lamb 
shift. 



fi-iTr{^ |(7,i3:|eff • E m\A\,^^)f'\ 

7 

pvr — i — -1 

L I K.% J 



7 

X PV[ ^ 



(7) 



In this formula {I7)} is the set of intermediate excited 
EM states of energy hio^ 



AB 



is the transition fre- 
quency of the ith atom, the script 5* stands for the time- 
independent operators in the Schrodinger representation 
and the script k can take the values 0, m, avg correspond- 
ing to an only dipole in free space, to mth configuration 
of fixed dipoles and to a statistical ensemble of dipole 
configurations respectively. Using an appropriate iden- 



tity for PV[ ^-T , it is shown in Ch.IV of 01 and 

[28] that the above expression can be written in terms of 
quadratic correlators as. 



equal point two-time quadratic correlation functions of 
the EM field and dipole operators in the corresponding 
EM vacuum and atomic ground state respectively [sol ]. 
In turn, the symmetric correlation function relates to 
the imaginary part of the linear response function via 
the fluctuation-dissipation theorem (FDT), 



dt exp [iujt]iQ{t) 



X {n,\[E{R\o),^{R\^t)]\n^)}, (10) 

(^j,|efi(w) ® ef-(a;)|Aj,) = -Tr-i^j J dt exp[iujt]iQ{t) 

x(A:,|[ef,(0),ef.M)]K)}, (11) 

and the whole problem reduces to computing the local 
linear response functions. That is, the polarizability of 
the dipole in its ground state and the Green function of 
the electric field, which can be calculated classically. In 
the following we elaborate on the physical meaning of the 
response functions. 



1. An only dipole in free space 

The equal point two-time commutator of the dipole 
operator for a two-level atom, say the ith one, in its bare 
ground state, \Al), prior to the coupling to the radiation 
reservoir and at zero temperature reads, 

I {eoh)-'eit){Al\[er'M,er'.i~mAl) 



^^{eoh)-'Q{t){At\ev,\Bl){Bl\ev,\Al) 
X [exp (iuJAB^) - exp (-iw^^)]. 



(12) 



The time- variable Fourier transform of the above equa- 
tion is termed the bare dipole polarizability, a{,(a;) -for 
brevity we have omitted the location dependence, 

a,{u;) EE 2(6oft)-^(^^|eFf |B^)(B^|eff |^^)a.$(y(^^^)2 - 



\ — 1 b,i 



= 2{eoh) uj^g 



b,i 



^ab) Ai^M, anyi, (13) 

- s 

(Aljefj jSj) is the unique dipole-transition 



cLSh _ 



dwK 



where ft = {AiL..^ ^j^^ 
^ -,. ^. matrix elemeno 6h,{lli) is the response function of the 

{Q.n\'Ei{R , w) (g) E (i? !'^)I^K)^jpQjg total monochromatic field acting upon i?' 



dt exp Mze(t) (Aj, I [ef, (0), eU~i)] } (8) 



r{Rr)^eoab{u^)-E'^,^{R'). 



(14) 



(A*Jer,(a.)®ef,(^)]A',) 



dt exp [iLot]iQ{t) {n^ \ [E(i?'' , 0), E(i?' , -t)] j(;9) 



where Eq.® is the energy shift associated to the polar- 
ization of the dipole due to the vacuum field fluctuations 
while Eq. (|77| is the energy shift due to the back-reaction 
of the polarized EM vacuum on the dipole. Therefore, 
the problem of computing f ^f* reduces to calculate the 



It is implicit that in order to get 'bare' atomic bound states some 
EM modes have already been integrated out in the definition of 

b.i 

The nomenclature used in this article varies slightly w.r.t. that 
in [23| |. In particular, o;, was denoted there as a' . Also, what 
was referred to as self-polarization field in is referred in here 
as radiation-reaction field. In addition, instead of referring to its 
propagator as polarization propagator, we name it here EM field 
propagator in order to avoid confusion. 
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For isotropic dipoles with a single oscillator -we can drop 
the dipole index, 



3eom; 



(15) 



where is the bare electron mass. However, since the 
effects of the EM vacuum fluctuations and radiation- 
reaction have not yet been incorporated, it does not sat- 
isfy the optical theorem. 

On the other hand, the time-variable Fourier trans- 
form of the electric field commutator in the EM vac- 
uum of free space, zeoft-ie(<)(rJo| [E(f, 0), E(r=', -i)]|f^o), 
is (w/c)^ times the Green function of Maxwell's equation 
in free space, G'*^°^(r — 



the atomic states to \Aq),\Bq) and, correspondingly, the 
single particle polarizability. Following a simple renor- 
malization procedure -cf. [Ij, [s^, Ull we obtain, 

d{uj)=do{l + ik^TT{ao-^[&°Hr,f-uj)]}^ \ (21) 

where the divergent real part of &^\f,'r;uj) has been 
incorporated in the renormalization of the free electron 
mass and the free space resonant frequency, nie , ui'^^ , so 



that ao 



Seam, 



Hereafter we will denote 



ujj^g simply by luq. a{u!) is the response function of a 
unique dipole in free space to an incident monochromatic 
probe field acting, say, upon , 



(22) 



^I - V X Vx 



;) = 5(3)(f-r=')I. (16) 



(5(0) (t?— r'-u}') is made of two contributions. These are, 
an electrostatic one, 



Gil.(^l^) 



PV 



r 1 

fc2 



(^-L)+k-'LS(^\r), (17) 



k = uj/c, and a radiative field. 



A kr 



—Aur 



„i kr 



(18) 



where the (5- function in Eq. (jl7|) must be intended in the 
sense of a distribution. The source tensor L takes ac- 
count of the geometry of the exclusion volume around 
each dipole-source. It satisfies TrjL} = 1 and, for a 
spherical volume, L = 1/31 [U. G(°)(r - f ; w) is the re- 
sponse function of the EM field to a unique point dipole 
oscillating at frequency uj at r' in free space, j^{'r'), 
E'^{r) — k'^tQ^G'^'^\r — f'-,Lij)-]J^{f'). In momentum space 
the radiative component is totally transverse w.r.t. the 
propagation direction while the electrostatic one is fully 
longitudinal. 



G(")(g;a;) = G^"' (g)P±(g) + G ((j)P|l (g) 



with 



— (f- 



1 



(19) 



(20) 



and Pi((7) = \ — q® P|i('z) — q® q being the trans- 
verse and longitudinal projectors respectively, with q the 
unitary vector in the direction of propagation. Hereafter 
and for the sake of brevity we will drop the explicit w 
and/or q dependence from the functional arguments un- 
less necessary. 

When the ith dipole, isolated, couples to radiation via 
the interaction Hamiltonian, the aforementioned radia- 
tive reaction of the EM field in free space renormalizes 



2. A specific dipole configuration 

Next, when all the dipoles are brought together, it was 
shown in [23j how a classical diagrammatic renormaliza- 
tion procedure leads to renormalized values for the single 
particle polarizability and the EM field propagator. In 
the first place, let us consider a fixed configuration of 
dipoles, m. The polarization propagator reads, 

N 

(23) 

where 7r,^(i?;,, i?^) is 



(24) 



= K-i% + fc2G(o)'(i?j„,i?^„)]-i. 

G(«)'(i?;„,i?^^) s G^°\Rl,,RiJ-^{G("HRla,RLm, is 
defined to take account of the regularization of the real 
divergence in &°\R\^,R\^). The inversion in Eq. (p4)) 
must be intended w.r.t. the particle indices, In par- 
ticular, the local term proportional to 5ij is defined as the 

ith renormalized polarizability, a^{uj). It is made of an 
infinite series of multiple-scattering diagrams which start 
and end at the ith dipole. Tri^{Rl^, R^^) is the (non-local) 
response function of the dielectric to a generic monochro- 
matic external electric field, 



N 



(25) 



The relation of 7r"j with the t-matrix defined in the 
so-called Coupled Dipole Method (CDM) [111 is tt^ = 

_l,-2Iw 

Next, let us compute the EM field propagator which 
yields the radiative corrections on al^{uj), gm(r, w). 
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It is given by an equation similar to that of Maxwell's in 
free space, 



;,,(f)i-vx vx 



(26) 

,=i,nS^'\^-RL) and e^/f) = 



with i-^^if) = 1 + ao Ej=i,Ar f'^^n?' - 
^mi^ ~ OioS^^^f— R\n). The source fixed at the position 
vector of the ith scatterer is removed from the permittiv- 
ity function on the l.h.s. of the equation. The radiation- 
reaction propagator is computed out of Qm making the 
source and the emitter coincide. It was calculated in [ll] 
in function of the polarization propagator, 



N 
3=0 



(27) 



and by consistency, 

~a'^ = ao[l + k^Tii^mil^^. Kn.^) ■ «o}]-\ (28) 

where again g„(i?j„, i?j„; w) = g„(i?j„, i?j„; w) - 
SR{G(o)(^J„,i?^;w)}. 

Applying reciprocity, gm(_R^, r*; w) yields the incident 
field which reaches a host dipole at J?^, whose source is 
a non-polarizable external monochromatic dipole at any 
point ^^^4, 



ext ■ 



(29) 



Note that since the host dipole at is polarizable, the 
total field at contains an additional contribution from 
the radiation-reaction (rr) field. As a result we have. 



From the above equation we can identify a new propa- 
gator which accounts for both the two-point, two-time 
commutator of the EM field in vacuum and the radiation 
reaction field fluctuations which dress up the renormal- 
ized polarizability of the emitter/source dipole, 

0m(^m,^;'^)= [l+fc^0m(^m,-Rm)-aoJ ' 0™ (i?^ , 7=') , 

(31) 



such that 



(32) 



g„j is the Green function of an equation like that of 
Eq. p6| but replacing ef^ j (r) there with the total dielec- 
tric function e"j(r). 



simply taking the statistical average on the precedent 
equations, 

N 

(33) 

or using diagrammatical techniques otherwise in momen- 
tum space -cf. [23( 1 ■ we end up with 

i ^ J dtd^rexp[i{u}t + q-r)]e{t) 

X {na.g\[E{6,0),E{r,-t)]\navg) 

= Qi- (9; uj){l-q'»q) + {q; uj)q®q, 



where 



x_L(g;^) 

pa 



G±{q;uj) 



t/||(q;w) = r G||(g;a;) = — 



X-L(g;^)/(P«) 

xii(9;w) 



pa 



pa fc2[H-X||(«';w)] 



(34) 



and Gj_j|(g;cLi) are the transverse and longitudinal com- 
ponents of the Dyson (bulk) propagators readily identi- 
fiable from the second equalities. The stochastic renor- 
malized polarizability, a, is defined later on. X-L,\\{q]^) 
are the transverse and longitudinal components of the 
electrical susceptibility. Diagrammatically, it is made of 
a series of one-particle-irreducible (IPI) scattering pro- 
cesses like those in Figslljb, c). The relationship of pro- 
portionality between the polarization and Dyson's propa- 
gator in q-space allows us to define the local field factors, 



pa 



such that II = Ci_ \\G 



^±,11(9) 

tively, Eg. dMl) can be written also as, 



.11^^.11 



Ahcr 



k'^pa 



G 



(0) 



(35) 



Equivalently, in terms of the stochastic polarization prop- 
agator. 



(36) 



N 



With Z^_2||:; (T) = / P^'[[d'RlS^'Hr-Rl^)S('Hr'~RL), 

1=1 

(37) 

we may write 



3. Random medium 



g^^^^{q;u:) = {p&r^Gf\{q;io)nl^^{q) 



(38) 



Finally, let us consider a statistical ensemble of config- 
urations of statistically equivalent isotropic dipoles. By 



For convenience we define a scalar radiative potential in 
terms of the trace of the radiation-reaction propagator, 



7 



(w) = + (j)'%uj), where, 



(27r)3 



-G 



(39) 



(40) 



For further convenience, we have distinguished the con- 
tribution of transverse from that of longitudinal modes 
in (j)'"^. The script sc stands for scattering since 4>'^'^{uj) is 
the part of (j) which involves multiple scattering processes. 
In function of 0, a and the renormalized stochastic po- 
larizability appearing in Eas. p4l38p read. 



a 
a 



ao(l 
q:o(1 



ao 
ao<; 



(41) 
(42) 



they all being isotropic, ao = ctol, a ~ al, a = al. 
When a single dipole is excited by a monochromatic 
external field, the emitted power is given by = 
^i^|£?^^jp9{c5;}, in agreement with the optical theorem. 
By parameterizing a and a as Lorentzian polarizabilities, 
formulas were found in 23] for the decay rate, F, and fre- 

2 

quency shifts. In particular, in free space, Fq = 377-^- 

As for the case of a specific dipole configuration it is 
possible to define a propagator which accounts for both 
the two-point, two-time commutator of the EM field in 
the stochastic vacuum and the radiation reaction field 
fiuctuations which dress up the renormalized stochastic 
polarizability. 



Q(r,r']u) = g{r,r';uj){l + a^c, 



(43) 



The total averaged electric field which acts upon a host 
dipole at f and whose source is an external monochro- 
matic dipole sited at r' reads. 



avg 



Sir, 



r 



H-ext- 



(44) 



For the specific computations in a random medium 
we will use a renormalization scheme similar to that of 
Felderhof and Cichocki [33]. It consists of two comple- 
mentary steps which refiect the mutual polarization of 
the dipoles and the EM reservoir. In the first one, the 
polarizabilities are renormalized by an infinite number 
of radiation-reaction cycles as outlined above. Diagram- 
matically, those processes amount to recurrent scattering 
terms in which the initial and final scatterers coincide. 
Note that those diagrams -cf. Figlljd)- may contain en- 
tangled intermediate recurrent scattering processes sig- 
naled by self-correlation functions. In the second step, 
the rest of IPI diagrams which are not accounted for in 
the renormalization of the polarizability are added up 
in the electrical susceptibility. Again, these diagrams 
may contain also entangled intermediate recurrent scat- 
tering processes -cf. Figsm&, c). Consistency between 



both steps is guaranteed by demanding that the scat- 
terers in all those diagrams be renormalized. Bearing 
in mind that unentangled recurrent scattering process 
are all accounted for in the definition of a, the electrical 
susceptibility components, ||(g;aj), adjust to cluster 
expansions of the form. 



E 

-l.m- 



X 



(n,m) 



n~n+m 



(45) 

The functions X^^'^^\q;uj) incorporate the spatial dis- 
persion due to the IPI spatial correlations within clusters 
of n scatterers in which all the self-correlation functions 
appear entangled. The index m stands for the total mul- 
tiplicity of entangled intermediate recurrent scattering 
events. The same kind of decomposition is applicable to 
e?_L,|| and (/.^^ -cf. FigUl 



— k pa 



(a) 




• • 



• • 



-• •- 



(d) 



FIG. 1. (a) Diagrammatic representation of Feynman's rules. 
Only two-point correlation functions are considered for sim- 
plicity. The self-correlation function, hseif{r) = p~^S^^\r), 
appears separated from the rest. (6), (c) Examples of IPI 
diagrams which amount to x'^' and x'^' respectively. The 
lowest order recurrent scattering diagram, x^'^'^\ is included 
in the series of x'^'- Examples of order p and or- 

der multiple-scattering diagrams which amount to C/'^' 
and (j'^-' respectively. The third and fourth diagrams in (d), 
which amount to Q^^'^\ are entangled recurrent scattering 
diagrams. In them, the self-correlation functions which affect 
the first and second scatterers cannot be factored out. 
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III. THE LAMB SHIFT AND THE TOTAL 
VACUUM ENERGY 

A. Lamb shift and Lamb energy in free space 

We are now in conditions to compute all the physical 
quantities. We start with the Lamb shift in free space. 
Applying the FDT on Egs- dSTfTl) with k = and substi- 
tuting there Eqs. p9l4ip we obtain the free-space Lamb- 

shiflU [il, 



was obtained by Buhmann et al. [ij, ll5| using a fully 
quantum-mechanical formalism. Making the identifica- 
tion 



AC 



'LSh,i 



= V 



N 



(46) 



From Eq. ipS)) we define the density of states contribut- 
ing to the free-space shift, AA(f'^''(w) = ^{^ao^^"^! + 

Next, using the Feynman-Pauli theorem -cf. pp. 295- 
297 of [33, the free-space vacuum energy is the energy 
gained by the system atoms-reservoir as the interaction 
Hamiltonian is turned on adiabatically neglecting dipole 
mutual couplings. This is parameterized by varying the 
coupling constant squared from zero to its actual value 
e^. Since ao is quadratic in e, we can write, 



(49) 

with the trace again intended over spatial indices only, 
the total vacuum energy density is 

7^1 JO Jo 

?{Tr[ln (% + k^ao ■ G(")'(i?j„, Ri^)^^ 



2ttV 



where the trace operation applies both over dipole in- 
dices («,j) and spatial tensor components. This result 
was firstly obtained by Renne (37j making the sum over 
normal modes and then by Agarwal [sHi using the FDT 
in an approach very similar to ours, and more recently by 
Emig et al. [38] using a variant of the Schwinger's source 
theory. Its expansion in multiple-scattering interactions 
is, 



n — 



hduj 



27r 



ph 

'2^ 



duj 3{ln [a/aolli 



dwS{ln[l-Hao(/>]} 
(47) 



result that was firstly obtained by Agarwal [36]. The 
free-space energies are additive and so we will refer to 
also as free-space Lamb energy density, Fq . 



Lamb shift and total vacuum energy of a fixed 
configuration 



Let us consider the mth configuration of dipoles. Fol- 
lowing analogous steps to those for E^^'^'' and J^g', we 
obtain instead, 



= F 



r°° ( ('„i'i"-i-i 

d.s{E^Tr{[. 

G(°)"(i?j„,i?^J]"}} 

/ Auj'il — Tr{(fc2a) 



27rV 



27rV 



(51) 



G(0)(i?:„,i?^„)}-f-Tr{(fc2«)3 ^ 



h 

2^ 

h 

2^ 



where G^"'" (i?:„, R^) = G(°) (i?L, ^^J - iK.. RDk 

In order to get a deeper insight into the physical inter- 
pretation of the above formula, we can write it in three 
equivalent manners. 



d^3{Tr[a„-0„,(i?:„,i?JJ]} (48) 
y{Tr[G(o)'(i?j„,i?^J • <(ff„,i?j„)]}. 



N 



2ttV 

h 



da;3i Tr 



dwSi Tr 



where the trace operation applies over spatial tensor in- 
dices only. The first expression after the equality symbol 



2ttV 

[G(")]-i(i?;„,i?^„J]]} 



lnK(i?^„,i?j„)/ao] 

In [0„«(^™,-Rm) 



}52) 



(53) 



2ttV 



du'^i Tr 



^ We emphasize that in the computation of any local physical ob- 
servable which is a combination of dipole and EM field oper- 
ators the inherent divergence in X{G^-^\RI^, W^; uj)} must be 
intended as regularized in both the free electron mass and the 
free-space resonant frequency of the oscillator. Hence, the use of 
regularized "primed" field propagators. 



In [a^ • 0„(i?*„,i?^) 
o-i.[G(°)]-i(i?^,i?^J]]}. (54) 



which shows that can be expressed in function of 
atomic d.o.f. only [Eq.dSH)], EM d.o.f. only [Eq.dSS])] or 
as a combination of both [Eg. ([5^ ]. In particular, the 
expression in Ea. ((55| in terms of the source EM field 
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propagator resembles Schwinger's approach [l7|. In ei- 
ther case, the normahzation by the energy of free-space 
fluctuations amounts to the substraction of the zero-point 
EM energy and bare atomic bonding energy. 



Average Lamb shift and total vacuum energy of 
a random medium 



Finally, let us compute the average vacuum energy 
of a random medium, J^^^ 0. In principle, since some 
stochastic observables can be calculated out of the op- 
tical response functions of the dielectric a{u}), fl", 
G'^, , it may seem possible to give a closed formula for 
•^avg which depends only on the electrical susceptibility 
and the renormalized polarizability. The computation is 
however far more complicated since the ensemble aver- 
age over involves highly non-linear terms in those 
functions. This can be expressed as 



Tr 



^Tr ln[(i 



In the first place, the ensemble average of Eq. 
the average Lamb shift, 



(55) 
is 



vLSh 
■'avg 



h 
h 
h 

2^ 



dwfc^cj 



(27r)3 
(27r)3" 



2Gi°'nj 



2^, 



doj S{q!0}, 



(56) 



(57) 



where we can identify the average density of states, 
TVf^^/i ^ pc>{|.0}. As noted by Bullough in [H, the 
integration of the equation analogous to that of Eq. ((50)) 
for J"^ 



= 

^ avg 



hduj 



^"O ^rLSh 

an 



(58) 



needs of the knowledge of the functional dependence of 
X on the polarizability ao. Thus, 
in closed form as function of x using Eq 



^otig'' can be expressed 



cLSh _ 
^avg 



h 

2ttP 



dwfc^S 



(27r)3 



From Ea. (j59l) we define transverse and longitudinal den- 
sity of states per unit of momentum volume, A/'iig''|i.||, 

k.ll =S{(2)^fc2xI||(9)Gf||[l + fc2Gf||Xi,l|]-^}, 

(60) 

where the prefactor (2)_l applies only to transverse 
modes, so that 



cLSh _ 
^avg 



hduj 



(2^)3[^-/U+Ae//l||]/p (61) 



However, except for some specific models, T^yg can- 
not be given in closed form. A way to go around this 
problem is to explode the cluster decomposition outlined 
previously for x, G and </). Firstly, by taking the en- 
semble average over the many-body expansion of Eq. (j5ip 
we can obtain a series for J-^yg in multiple scattering 
terms. Note that this is not yet a series in p since, 
at a given order n of Eq. (j5ip . there are processes with 
2,3,..n different indices. When performing upon them 
the ensemble average, they contribute with terms of or- 
der pa4>^^\..pa(j)^"'~^^ respectively. Equal index terms in 
Eq. ([5T|) involve a self-correlation functions 6'^'^\Ra — Rb) 
when performing the ensemble average. 

Further, some realistic approximations can be carried 
out in order to obtain an expansion in powers of p. In the 
first place, radiative-reaction corrections in the renormal- 
ized polarizabilities which enter any diagram in can 
be disregarded in good approximation. Upon integra- 
tion in w both in Eq. (|57| and in the average of Eq. ([5T|) . 
they yield terms of order ^ fco^'e ^ Fo/wo ^ 1 smaller 
than the computations involving bare polarizabilities, 
being the electron radius and Fq being the free-space de- 
cay rate of an oscillator. Second, entangled intermediate 
recurrent scattering events can be ignored since they pro- 
vide terms of orders [(/coO'^To/wo]^™ <C 1, m > 2, s = 
1, 2, 3, .., n, smaller than the non-recurrent, non-retarded 
near field terms and terms of orders [(fco^)~^Fo/a;o]'"" ^ 
1, m > 2, s — l,2,3,..,ri, smaller than the non- recurrent, 
retarded radiative term at any given order n in p, with ^ 
being the typical spatial correlation length, s the number 
of scatterers repeated and m the recurrent multiplicity - 
cf. computation of (/i^^^ in SecfVland Appendix VK\ With 
these conditions provided, we can simplify Eq. (j57p and 
-ithe average of Eq. ((5T|) as 



+ X||G[°)[1 



fc^G[°)x||] 



LSh 



(59) 



^0 



LSh 



cLSh 
^0 



?LSh\ 



h 

2^ 



LSh 



{2n) 



E 



,(m,0) 



In principle, it could be possible to include the randomness of 
the dipole configurations in the actual (not averaged) internal 
dynamics of the dipoles as considering the computation of the 
actual (not averaged) total vacuum energy. For this to be the 
case the correlation time of the density fluctuations should be 
much less than the typical dynamical scale, ~ 2it/ujo- However, 
under these circumstances it is likely that the dynamics of the 
external d.o.f. be fast, which would induce additional complica- 
tions like Doppler effects. 



avg 



■^avg IctQj 

h 

2^ 



} 

(62) 



■J n 



oo 
m— 1 



.0) 



} 



avg \oco ■> 



(63) 



where the scattering Lamb shift, 

,(m,0) 



explicitly and 



^^'^ has been isolated 
(uj) stands for the m-body term of 
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the expansion of (j) in which no entangled intermediate 
recurrent scattering events appear (norec.) and polariz- 
abilites are taken as bare. 

Lastly, we can write Eqs. (l62l63p in function of the clus- 
ter expansion for x- This way, we obtain a formula for 
J-avg\aa which depends explicitly on (the cluster expan- 
sion of) X, order by order in p. As for Ea. (l59p we use 
Eg. (1341) and the precedent approximation to write. 



they are model dependent. Nonetheless, they are useful 
for estimating orders of magnitude. In particular, that 
in [23] corresponds to the so called quasi- crystalline ap- 
proximation (qc). According to this approximation the 
only relevant correlation function is the two-body one, 
/i(r), and self-correlations are ignored. In this approxi- 
mation, the series of x(9) becomes geometrical and the 
only quantity to be computed is x'"^''^' (9)1 



2Xj 



norec. y^(0) 



(27r)3 

(-1)"+! 



[1 + ed^^xTZ-r x^^A) = I dVe--^MOTr{G(°)(f) • P,(^)} 



2y^(0) ^^norec. i — 1 

OLQ \ 



(64) „ r. ^^^^ 

with p =±, II . The geometrical series in Figi2]is 



(27r)3 ^ 

^ ' m=l 
2 ^^norec. ^(^) 



2[k'xr 



\,ao 



(65) 



(69) 



where x"°p„|^ are the transverse and longitudinal compo- 
nents of the susceptibility with bare polarizabilities and 
no entangled recurrent scattering processes. Using the 
cluster expansion of Eq. (|45l) we obtain. 



cLSh\ _ cLSh 
^avg \ao ^ '-'0 



— / da;3{ 
27rp ' ^ 



(27r)= 



-(fc^x^:°!Gf)^ + (fc^x^:°;Gf)3 

+ fc^Gfx£°)-2xt°US!,(fc^Gf)^ + ..] 

+ fc^G[,^[,^;,°^2x[,::^x£^fc^G;,°^)^ + ..]}, (66) 



h 

2^ 



(27r)3 

-(fc^xK!Gf)^)+i((fc\(-)Gf)3 



rl 



2r^ii,"o^ii ^'^ ^ 

2, (1.0)^(0)n3 



,2,12,0)^(0) 



2,/1.0)^(0)n2 



2^(0)^/3,0) 

ao 



2x[,::'xi,-i(fc^G[°))^) +..]}, 



(67) 



which in turn yields an expansion in powers of pe^ like 
that of Mclachlan [39]. It is plain that the problem in 
integrating Eq. ((55)) has been shifted to that of knowing 
X_L.||(<z;w) as a power series of p or, otherwise, knowing 
all order spatial correlation functions. 



IV. APPROXIMATIONS TO THE TOTAL 
VACUUM ENERGY OF A RANDOM MEDIUM 

A. The quasi-crystalline approximation 



In previous works -cf . [8|, [231, |2J] , it has been reported 
closed formulas for J-^n as functions of x±,||- However, 



and so, 



V(2,0) "\ 

yia "0 — 



X,, 



0--kpa 

0+6 — 6+6 — 6 — 0+6 — 6 — 6 — 0+.. 



-6—6—6 



-6—6- 



FIG. 2. Diagrammatic representation of Ea. (|69 

(9, ^) - Gj; [1 + pao(fc^G5| - x^;L/p^«^)]-^ 



(70) 



and 



nr 



qc 



2tt 
X In 

h 

~ 2^ 



qc y2gqc 



(27r)3 

„JGf]-[Gf] 



} (71) 



(27r)3 

X '2ln[l + pao{k'Gf-xfZ/p"^o)] 

+ ln[l + pao(fc^G}°'-x£°Vp^«g)]]}, (72) 

where by ignoring radiation-reaction it holds that = 
^[[^l|. It is because of the aforementioned analogy be- 
tween the formulas for a specific configuration of dipoles 
-cf. a cubic lattice @- and those of the quasicrystalline 
approximation, that the inequality expressed in Eq. ((55)) 
turns into an equality in this case. It is also possible to 
write Eq.dni) as 





2n 



(27r)3 
d^q 



21n[xrGl1+ln[xfGf 



(27r) 



2 In 



[G 



(0)] 



In 



[g[°^]] ^74) 
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Making the identification J ^-^^^ = p in Eq. ([75l) we 
read that the bare atomic bonding energy and the EM 
zero-point energy enter as substraction terms. This jus- 
tifies the interpretation of J"^ in Q. That is, J"^ takes 
account of the zero-point energy of bare EM modes in 
Eg. ([74|l ■ and the atomic bonding energy of Eq. ([75l) and 
substitutes them with the binding energy of the coupled 
system of Eq. ([73ll . 

Further, by considering the continuum limit of the 
effective dielectric constant, ^ 1, we will see that 
^ 1) can be written as a function of the refractive 
index. We will discuss the accuracy of this approxima- 
tion. 

Although the expression obtained in [l^ for J-'Yvg co- 
incides with that of Eq. (iTlT) , the derivation there was er- 
roneous and hence its validity restricts to the quasicrys- 
talline approximation. The reasoning followed in [23| was 
that the average energy could be obtained by extending 
appropriately Schwinger's approach on an effective con- 
tinuum medium [19] to a molecular dielectric. Thus, the 
ste ps f ollowed by Schwinger et al. in [i^ and Schwinger 
in [iSl [T9I were mimicked in [2^ but for the fact that 
the effective bulk propagator was replaced by the aver- 
age source field propagator of a molecular dielectric. The 
reason for doing this was the constatation in that the 
effective bulk propagator used in [lEIiSllSl is not the one 
which enters the formula for the Lamb shift in Eq . ([56| . 
Accidentally the approach in 23^ turns the inequality of 
Ea. (|55p into an equality, which is generally incorrect. 



B. The Schwinger bulk energy of an effective 
medium 

The reason why the ap proach of Schwinger et at in 
[ioj and Schwinger in [3l does not yield the total bind- 
ing energy of a molecular dielectric is that the variation 
of the energy shift considered there is not appropriate 
to this aim 0. The EM effective action of Schwinger et 
al. [40| is that of a piecewise homogeneous medium of 
effective electrical susceptibility Xeff (^) with interaction 
Hamiltonian, 

H:J^/ = - j d\P,ff (f, t) ■ Eij f (f , t) . (76) 



This equivalence can be proved passing to the continuum from a 
lattice of dipoles whose cell volume is ■ It acts as a regulator. 
Note that this does not imply by any means that the generic 
Schwinger's source theory be inapplicable to this aim as shown 
explicitly for the case of a specific configuration of dipoles, 
Eq.l|53} -cf. O. 



The corresponding energy shift is 

SSch. = -5R{ y dVd {nSch.\Peff{r,t) ■ E^^^(f,t)|175ch.)} 

(77) 

where E-^f f is the transverse effective field generated by 

an effective polarization density source Pef / through the 
relation, 

E^f^{f-u;) = kh^' J dVGl^/(f,r';c.).Pr/^(r-). (78) 

The response function of the source field is the trans- 
verse effective bulk (Dyson) propagator, C^f^ [q-.uj) = 
{ceffk'^ — q^)~^ ■ Adapting the previous nomenclature, 
G^/^{q;Lu) characterizes the EM field fluctuations of the 
'Schwinger vacuum' (Sch.). By var yin g the action of the 
effective medium, Schwinger et al. [40| inferred the iden- 
tification of the effective polarization source quadratic 
fluctuations, 

J dtexp[iujt]i{eohy^e{t) 

X {nsch.\'^eff{r,0)(^'^eff{r',t)\nsch.) = SXeffS'^^'Hr 

Inserting the above expression together with Eq. ([751) into 
Eq.dZZl) we get [il, 

x[l + k'GTx'"„]-'}- (79) 

The functional integration immediately yields 

^^...^^f d.3(/^l„[R"P|G™l-],, 

(80) 

which contains only the fluctuations of the transverse ef- 
fective bulk propagator. 

1. Relation with the microscopic approach 

In order to understand the link with the micro- 
scopic computation we first notice the similarity between 
Eq. (j80|) and Eq. ([59| for the microscopic computation of 
^avg^ ■ Ignoring the fact that Schwinger neglected the 
contribution of longitudinal modes, Eq. (j80[) can be writ- 
ten also as 

where N^^f^'f^ is the average density of states re- 
stricted to transverse modes in the continuum. This 
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implies that, effectively, in the approach of Schwinger 
et al. the microscopical interaction Hamiltonian of 
Eq.(IT]) with dipolc moment operator ef and coupling 
constant e is substituted with the effective Hamiltonian 
of Eq. (l76|) . which couples the transverse effective EM 

field to an effective polarization density operator Pg// 
with some effective coupling constant proportional to 

1/2 

Xeff- Because Xeff is made of the integration of dipole- 
dipole interactions, the energy of those interactions is 
disregarded in the approach of Schwinger et al. Q- This 
means that Schwinger et al. approach is not appropriate 
to study phenomena driven by the variation of the 
energies of internal d.o.f. This is the case of Lamb shifts 
van der Waals forces and phase transitions in liquid 
crystals [4l[ and colloidal suspensions 42 1 in which the 
relevant variations in energy are those w.r.t. that of the 
molecular constituents infinitely far apart. However, it is 
appropriate for the treatment of problems in which the 
internal properties of the objects are integrated in the 
definition of their permittivities and remain unaffected 
during the phenomena under study. This is the case 
of the Casimir forces between macroscopic dielectrics 
separated by macroscopic distances [40] in which the 
relevant variations in energy are those w.r.t. that of the 
macroscopic objects infinitely far apart. Nonetheless, 
the reason why Sq^'^ can be computed out of J^Jc/i. 01 
is that, at leading order in p, x(r) ~ paolS^"^") (r) and so 



fV 



J'lh. at Oip). 



Sch 

The original suggestion of Schwinger to explain 
sonoluminescence represents a sort of intermediate 
problem [1^. Basically, when a bubble embedded in 
water collapses, molecules of water fill in the void and he 
conjectured that the light emitted in this process carries 
the excess of EM vacuum energy stored in the void 
w.r.t. the energy of the homogenous aqueous medium. 
Evidently the molecules of water which fill the void 
form part of the initial aqueous medium surrounding 
the bubble. Therefore, leaving aside the validity of 
further dynamical approximations - cf. [43| . the total 
volume filled by the water is greater after the collapse, 
the density of water is less and so the dielectric constant 
must vary, at least locally, in the region within and 
around the primordial bubble. As a consequence, it 
is reasonable to think that, if the time scale of light 
emission is less than the typical homogenization time, 
the molecules of water within the volume originally 
occupied by the bubble will have different spatial dispo- 
sition to that of the homogenous surrounding medium. 
Therefore, their internal energy will be different to that 
of the homogeneous phase and the problem involves not 
only a variation of energy between macroscopic objects 



''' Analogously, in assuming a given coupling constant and a bare 
resonant frequency, the computation of the internal binding en- 
ergy of the atomic dipoles prior to the coupling to radiation is 
also disregarded in our approach. 



(aqueous media with and without a bubble) but also a 
difference in the internal energy of the molecules which 
were initially in a homogenous disposition and fill in the 
void afterwards in a different arrangement. On the con- 
trary, if the homogenization is reached prior to emission 
and, in good approximation, the permittivity of the filled 
bubble is identical to that of the surrounding medium, 
Schwinger's approach might be a good approximation. 



2. Beyond the continuum approximation 

The result of Schwinger in the continuum yields only 
the energy associated to radiative modes which de- 
pends solely on the effective refractive index, n(a;) ~ 
as 



-> Sch. 



3? 



and the Schwinger energy shift is 



(82) 



(83) 



In Sec lVBI we compare in detail the above results with 
the actual microscopical calculation. Schwinger argued 
that these results in the continuum could be ex- 
tended to incorporate spatial dispersion in the electrical 
susceptibility. In doing so we obtain an extended (ext) 
Schwinger transverse bulk energy, J^^^y^ ]ext ^ However, 



even in this case the discrepancy between the actual J^^ 
and ^Jc/i 1*^^* appears already at order p^ . For the sake 
of simplicity we restrict ourselves to the approximations 
used in Ea. (|571) and find. 



^ avgiOio — ^ avg idQ 



-pV \ext I 
Sch.lao + 



2tt 



d^g 
(2^ 



1 



.2,12,0)^(0) 



l,2M0)im ^ (1,0) (2,0) 2^(0)^2 



(84) 



where J-, 



V ill 



contains the longitudinal modes whose 
contribution includes radiative energy -cf. computation 
of -^o(p2) in SecES 

For the sake of completeness we mention that the ap- 
proach of Abrikosov, Dzyaloshinskii, Gorkov, Lifshitz 
and Piaevskii (ADGLP) [H, is based on a semi- 
phenomenological prescription for the vacuum energy 
which accounts for the restriction of the EM fluctuations 
to those which contain a strictly local polarization oper- 
ator, ^ Xe//I^^'^' (^ ~ ^)- Hence, its content is similar 
to that of Schwinger's in the continuum but for the fact 
that it includes the contribution of long- wavelength bulk 
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longitudinal modes 

^ADGLP = 



doj 
2^ 



(2^ 



In 



ee//(ee//A:2 -92)2 
.(85) 

where e^ff is the effective dielectric constant. A critical 
analysis of it has been carried out by BuUough in [i^l • 



V. THE LEADING ORDER LAMB SHIFT AND 
VACUUM ENERGY IN A RANDOM MEDIUM 

At leading order, the free-space Lamb shift and Lamb 
energy are additive. Also, at leading order in p, ~ 
pEq^^- The actual computation oi Tk , £ n ^^ has been 
carried out by a number of authors [Tu [211 [U, , 



r f /-"^ do; ^ uj^ 



(86) 
(87) 



dw 



(27r)2c3 



Upon substraction of the divergent free-electron self- 
energy [2l|, = ^j^TcT /dw w, the above integral 
presents a UV divergence which needs to be regularized. 
In the non-relativistic approximation, the 'natural choice' 
for the UV cut-off A is that of the Compton wavelength 
of the electron such that A^ — 2-kc\q — uieC^/h. For a 
single oscillator with = 



J-r. ~ — at In 



hujo 



2tt 



In 



(89) 
(90) 



where 



Bethe's result when expressing the atomic ener 
terms of the velocity of a bounded electron [23, Ea] El- 
Should have we computed Eq. ((86| and set the wavelength 
cut-off at the electron radius, AJ; = c/r^, we would have 
obtained the classical result of Bowling [13] by consid- 
ering the singularity around this value. However, this 
is inconsistent with the non-relativistic approximation. 
The need of a cut-off just reflects our lack of knowledge 
both of the internal structure of the dipoles and of the 
manner the EM field couples to the internal d.o.f. 



IS 



is the electron radius and ctf = 

Eq.dMl) equals 
in 



47T€om 

the fine-structure constant 



Also, a similar expression to that in Eq.| |89I | was obtained by 
Welton [T3II from the variation of the non-relativistic Compton 
scattering cross-section due to the position fluctuations of the 
electron. 



At higher orders the Lamb-shift is non-additive. At 
0{p) we have. 
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^a0«(-)}, (91) 



where the subscript a means that the renromalized por- 
larizabilities in (l>^^^{uj) must be replaced with free-space 
ones in order to keep the order p of ^ofp)- ^^^^ order, 
the corresponding 0- factors read from Fig[T]((i), 



Q-L. 



(27r)3 



-pa[Gj'||]^+xL'i,l|G5|/(pa) • (92) 



From the above equation we can write, in function of 
LFFs, 
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hp'^ 



duj 
2^' 



(27r)3 



xS (<?)/(/>«)' 



2eGf\q)[x':iiq)/{pa) 



(2) , 



(9)]]}- (93) 



This is to show how LFFs enter the Lamb shift at order 
p. Related to this fact, the authors of have computed 
the effect of LFFs on the van der Waals forces on an only 
dipole in an Onsager (real) cavity using Eq. (|9T|) . 

At leading order in ao, it is plain from Eqs. (|62l63p 
that J'o(p2)\ao — P^ofp') Iqo/2- From those equations the 
same simple relation holds at higher orders, J-Y 



p_cLSh I 
m''e)(p™-i)|ao- 



O(p'")l"o — 

This allows for an expansion of both 



J^avgUo and £av!^\aQ in TO-body terms of order (e^p)™. 
However, beyond the approximations used there, no sim- 
ple relation exists since any given order p™ contains 
higher powers of due to entangled recurrent scattering 
and additional radiative corrections. 



A. Computation in the hard-sphere model 



LSh 



Except in free space, the actual computation of Ej^^g 



and FY„n is model dependent. 



Nonetheless, 



generic 

results can be obtained within the simplest analytical 
model. Let h{r — ^) be the two-point correlation func- 
tion, ^ being the correlation length. Generally, h{r — ^) 
can be modeled by the addition of three terms, 

h{r - - /i'"'- (?- - + P^^f^^^' {r) + h'""^{r - i). (94) 

In this equation, h'^^'{r — ^) accounts for the exclusion 
volume around each dipole and its precise form depends 
on the interaction potential between pairs of scatterers. 
It tends to -1 for r < ^ and to zero for r > ^. Usual forms 
are those of a Lennard- Jones potential and a hard-sphere 
potential. h°'"^^{r — takes account of the overdensity 
of first neighbors around a given dipole. It might be 
relevant for high-ordered media. The three-dimensional 
delta function stands for the self-correlation, h^^^ and 
the self-correlation functions are inherent in any molecu- 
lar dielectric. 
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In order to make contact with previous approaches we 
further demand the existence of an effective medium for 
the frequency range of interest. That imphes C = A:^ ^ 1 
for u! ^ ujQ. For the sake of simphcity we wih neglect 
in first approximation both the overdensity and the self- 
correlation terms in h{r — We will show a posteriori 
that the latter is plainly justified while the former implies 
slight modifications in numerical prefactors. 

Without much loss of generality we will take a 
hard-sphere (hs) exclusion volume correlation function, 
h^^-^r — ^) = /i'"'(r — ^) = —Q{r — which derives 
from the potential U{R\ ..,^^) oo if \R' - < 
i ^ j, and otherwise. The numerical factors of the cal- 
culations which involve near field modes depend on the 
precise profile of h'^^\r — ^). On the contrary, radiative 
propagating modes are model-independent. 

The computation of in spatial coordinates is eas- 
ier than that of Eq. (|92|) in Fourier space. 



/,(1,0) 

Q.,llS 



= Tr 



d^ 



G("'(f)(- 



-fcVa)G(°)(r) 



X [1 - e(r 
pae 



27r 
27r 



■pa 



0\] 
1 

1 



6 ^ 



(95) 



C < 1. (96) 



The decomposition into transverse and longitudinal com- 
ponents is given in Appendix \K\ Inserting Eq. (j95p into 
Eg. (1911) and integrating in uj we obtain. 



cLSh \norec. 
■'0{p)\hs 



12eo ujo 



(Co^'-Co"') + ^(5/6-7£"ln[2Co]) 



(97) 

where j£ is the Euler constant and, as for the compu- 
tation of J-Q, the cut of the integrand at has been 
neglected. For simplicity, the integral has been expanded 
in powers of Co = fcoC ^ 1 up to order zero, b eing the 
leading order term that of London's potential [49[. In 
contrast to J^q, the oscillating factor e^*"" serves a natu- 
ral UV cut-off at A ~ c/2^B 

In the last equations the superscript (1,0) signals that 
recurrent scattering terms have not been included. As 
advanced in Sec lIII C[ the condition Co » Fo/cjo ~ fco'^e, 
implicit also in Ea. (|97l) . suffices to guarantee that recur- 
rent scattering can be neglected in good approximation. 
In Ea. (jB3[) of Appendix IB] we give the expression for the 
complete series of recurrent scattering diagrams which 
amount to (f'a'^ . Its expansion in powers of a and further 
integration in uj yields a series of the form, 

.To 
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/™^'(Co)(roM)^ 



Wo 



® Making use of the property du}Q{(f>{uj)} = du}<f>{iu)) 
[I2l |36|1 . all the integrals in uj are performed with the change of 
variables u = —iuj. Since our computation refers to the ground 
state, resonant terms are absent. 



The functions fm^^iCo) contain both negative and pos- 
itive powers of Co together with terms proportional to 
In [(2(m + l)Co]. Each order m presents a wavelength 
cut-off at 47r(m + 1)^. For m = 0, fo^'^iCa) is readily 
identifiable from Eq. ((97l) . The leading order term of a 
generic function fm^^iCa) is of the order of ^~3(2m-hi) 
Therefore, [f^S\{r,/u;o)'^"^+'^/[f^f'^{ro/uJo)'"'] ^ 
/ ro/i^o N2 convergence of the series is guaranteed 

''0 

by the aforementioned inequality, Co ^ To/'^o- At the 
same time, this makes the neglect of recurrent scattering 
terms a good approximation. 

Because each order in recurrent scattering carries an 
additional factor Oq, the integration of Sofpjlhs yields for 
non-recurrent and recurrent terms respectively. 



J", 



\7iorec. 
\hs 



'PP 



24eo cjQ 
-ln[2Co]) 



2 r r 



dip^)\'^s^P~e,^p~-^Y.— 

ni—O 



-rV I 2^-1, ,2 



/™"''(Co)(roM 



is 



(99) 



,2771 



(100) 

The numerical constants in Eq. ([77)) depend on the spe- 
cific profile of h^^-. Nonetheless, the order and the sign 
of the terms are generic. Hence, the script hs can be 
dropped from Eas. (|98ll00p since those expansions are 
not constrained to any particular model. As a matter 
of fact, the addition to Ea. (|n5|) of the terms in Ea. (|A9|) 
which account for the overdensity two-point correlation 
function, h'""^{r) = C^S'-^^r - 0, just modi fies sli ghtly 
the numerical prefactors of the terms in Eqs. (|97l99|) . but 
neither their order nor their sign. 

It is worth mentioning that, although at a given or- 
der the leading contribution comes from non-recurrent 
terms, there are recurrent terms in orders p*, < s < n, 
which are of order (pC^)"*-""*' > 1 greater than the lead- 
ing order non-recurrent term of the order p" . 



B. Discrepancy between J-Xvg\rad and T^^h. 

We investigate to which extent the Schwinger vacuum 
energy of an effective medium is a good approximation 
to the actual vacuum energy of radiative modes. We 
will compare their relation up to order p^ . For the rea- 
sons given in the previous section recurrent scattering is 
negligible and the effective susceptibility is well approx- 
imated by that of a Maxwell-Garnett (MG) dielectric. 
An MG dielectric is characterized by the fact that the 
only relevant correlation function is that of an exclusion 
volume, h'^^-{r ~ £). In the electrostatic-long-wavelength 
limit it is proven that the quasicrystalline approximation 
is exact [501 [^^e FigO] and the effective susceptibility is 
independent of the precise form of h'^^\r — ^). In par- 
ticular, we can use the results of the hard-sphere model 
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without loss of generality. According to this, the elec- 
trical susceptibility is the sum of a geometrical series of 
ratio kmqI pci — potji, in which only the bare longitu- 
dinal propagator and free-space polarizabilities enter. It 
yields, 



Xmg 



pa 



pa 



(101) 



1-X^Wpa l-pa/3' 
and from here and neglecting further renormalization 



0=-k^pa 



and so J-g^f^ P*-'' — -7^Q(p2) Irad- However this equal- 
ity is accidental. To see this, it is necessary to express 
J-^yg and J-^^j^ in comparable terms. To this aim we 
use the expression of Eq. ([63|) for J^^g but with a instead 
of ao there, and the extended version of J^sch which 
incorporates spatial dispersion [31, -^g^^^ 1*^^*. Inserting 
Eq. ((35)) into Eq. ([80)) for J-g^i^ but with spatial dispersion 
we have. 



■j^V \ext 
■f Sch.\ 



d^g 
(2^ 



In 



\-pa J2 '-'^^"^ 



Further, we expand the latter equation up to order 



(104) 

p'B 



FIG. 3. Diagrammatic representation of the IPX processes 
which amount to the MG susceptibility. Only the exclusion 
volume two-point correlation function, h'^^'{r — ^), is relevant 
and self-correlations are disregarded. In the electrostatic- 
long-wavelength limit, k^, 0, the quasicrystalline ap- 
proximation is exact. 

on a for the reasons given above, the effective refractive 
index is n ~ 1 + pa/2 + (pa)'^ /24+ ... It has been already 
seen that J-^g^i^ contains the free space Lamb energy. At 
0{p'^), by inserting the series of n in Ea. (l82|) we get, 



Sch.l 



7h 
'48^ 



5ft 



did t 



■■{pa)'}. (102) 



Next we turn to the hard-sphere model for the micro- 
scopical computation. The C- independent terms of (j^a hl 
in Eg. (1961) amount to the energy of radiative propagat- 
ing modes (rad) within J-Qi^p2f\rad upon integration in 

oj. We show in Appendix [K\ that, out of them, —i^^pa 

comes from the bulk transverse propagator in 2^'^'°'' F"l 

[Eqs. (|Al|A3IA7p ]. The remaining -i^^ comes from 
the bare transverse propagator in the susceptibility func- 



tion of 



[Eqs. (IA3IA8|) ]. Hence, we can write 



7h 

'487r= 



3? 



dw t 



■'{paf}, (103) 



V lext 
Sch.l 



'■{'ia^iVp- 



du) 
2^ 



d^g 



[g|l°¥ (105) 



2a^k^ 



-a^k^ 



{2tt 



^GfV-^^/p + Oia^))]^6) 

On the other hand, for the evaluation of the radia- 
tive modes of J^avg use the expressions for the MG 
factors and effective transverse propagator restricted to 
radiative modes ^23^, 



gf^((z) = /:LLGl{^^(Q), <PmG 



2tt 



Cl^n, (107) 
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LL 



n, 0A/GII = (f'MG - 2(j)MG±, 

(108) 

[{xmg + l)k^ - 



and Gl{^^(g) 



with Cll = 

From Ea. (ll05p we have that the first term there con- 
tains only the transverse modes in J^Q(^p2y Its restriction 

to radiative modes amounts to — 2l§^^{/o°° '^^ k^{pa)'^}. 
On the other hand the second term of Eq. ()105p amounts 
to i^^if^ d"^ k^{pa)'^}. Adding up the last two quan- 
tities we obtain Eq. (|102l) which equals J'^^p2^\rad, even 
though no longitudinal terms enter the Schwinger result. 
The reason for this equivalence is accidental, since the 
0{p'^) radiative term of the difference J^Yvgl^ ^ -^Sc/i.T"^* 



A comment is in order here to amend some erroneous interpre- 
tations in |23|I . In the first place the distinction between co- 
herent and incoherent radiation carried out in Sec. IV of [23|l 
is erroneous. While Ea.[23jf54') is correct, its equivalence with 
Eqs. [2^ (42+44,46,48,50) is not. Consequently, the conclusion 
that only one LFF appears in the expression of the coherent ra- 
diation is incorrect. The correct calculations will be published 
somewhere else |5H] . 



The expansion of the logarithms in Eq.l(T04} yields a series of q- 
integrals whose integrands are products of powers of polarization 
propagators of the form [S'^l]". Each factor is an m- 

scattering loop and the sum of the exponents, s, is the loop order 
-eg. G^'ii [S^'iiJ^s'^'^'ii is a four-loop term made of one free-space, 
one single-scattering and two double-scattering loops. 
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in Eg. dM)) is equivalent, by reciprocity '23'l, to the longi- 
tudinal term in /"^jp2)| rad but with opposite sign. This 
relation is not model dependent. Note also that by incor- 
porating spatial dispersion in J^g^f^ |'^'^*, the term 20^'"^ 
contains only one half of the non-propagating term in 
Eg. dMl) proportional to (q^. The rest of near field terms 

^0{p^) are disregarded in Tsch.loip^)- 

Next we show that the accidental coincidence of 



J". 



V text 
Sch.\o{p^) 



and J-, 



breaks down at 0{p^). Ei- 
ther by expanding the integrand of Eq. ([82l) or by direct 
integration of the 0{p^) terms in Eq. (|106l) we get, 



V \ext 



Sch 



I7h 



duj k 



■■[paf], (109) 



2887r2 .j^ 

while the insertion of Ea. (ll07p into Eq. ((63|) yields. 



17ft 
1447r2 



dw k 



\paf} 



= 2 J", 



Sch.\o{p^) 

(110) 



VI. THE LAMB SHIFT FROM SCHWINGER'S 
EFFECTIVE MEDIUM VACUUM ENERGY 

Schaden, Spruch and Zhou (SSZ) '22'] have computed 
the Lamb shift at leading order in p from the Schwinger 
energy of an effective medium in Eq. ([5^ . i.e., the free- 
space Lamb shift. For the reasons explained in Sec lIV Bl 
Schwinger's approach yields the correct result at lea ding 
order for Xeff — po^o- Nonetheless, the authors of |22| 
have applied a more general procedure consisting on com- 
puting the Lamb shift out of the variation of -Fg^i^ w.r.t. 
"small" variations of a background refractive index. That 
is, adopting our nomenclature, they have used the vari- 
ation of Eq. ([75)) together with the result of Ea. (|551) and 
Sx 



2nAn, 



^SSZ — P ^n-r Sch 



poo 

(Ill) 

For a dilute medium, n ~ 1, /S.Xeff — pcta, and 
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^ssz 



-h 



3?{ / dww^pao}, (112) 



which is nothing but the free-space Lamb shift. This 
is analogous to the com put ation carried out by Feyn- 
man. Power and Milonni [111 IMUl in the limit pV 1 
in which the medium consists of an only dipole in free 
space-V being the sample volume. More specifically, the 
authors of 22] have computed the difference of the free- 
space Lamb shift between two dielectric states, / and //, 
with Axeff = p{al' - ai). 

The computation of Milonni, Schaden and Spruch 
(MSS) in (3 combines elements of both [ill, [13, and 
[23 |. That is, while in Q the variation of the bulk en- 
ergy density is calculated between two different states of 
molecular dielectrics as in [l^l , the difference between the 



states is given by the difference on the polarizability of 
only one of the dipoles as in [ll], US IHI • Under the as- 
sumption that the dipole in question is randomly placed 
and the medium behaves as a continuum of refractive in- 
dex n ^ pa^'^^ /2, the authors take Axeff = p(ao^ ~Q^o) 
and, using Eq. ()llip . they get 
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/ dLUio^np{al,^ ~a'^)}. (113) 



Likewise, for the Lamb shift due to the presence/ absence 
of the background medium (BGM), Sa/ssIbgm ^ 
d^(a.Vci)(n - iKai' -ai)}. 
The authors of had already warned that the deriva- 
tion of the Lamb shift this way might need to be corrected 
by local field factors in high ordered systems. We have 
proved in Sec. IV Bl that this is indeed the case under 
any circumstance, no matter the degree of order. As it 
was mentioned there, spatial dispersion and longitudinal 
modes are to be added to Schwinger's formula in order 
to obtain the correct result. Even if the host medium 
can be treated as a continuum, the microscopical calcu- 
lation of the Lamb shift on the simplest geometries for 
the embedding of a dipole contains LFFs. At leading or- 
der, near field and other radiative factors enter the Lamb 
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For instance, for 



shift through the terms oi m '--Q(^py 
a small Onsager cavity (0ns.), those terms are the ones in 
Eq. dHU) with ^ being the radius of the cavity and kS, < 1 
[H, m, HI- Thus, instead of Sm^Ibcm ^'^ obtain at 
leading order in (n — 1) 
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(114) 



— 3{/dc.(n-l)(a^^-a5)[^ + ^]}. 



VII. CONTINUUM APPROACH TO THE 
BINDING ENERGY OF AN EFFECTIVE 
MEDIUM 

In this Section we investigate the possibility of quan- 
tifying, at least partially, the binding energy through 
optical observations. This is motivated by the conjec- 
ture raised in the introduction on the correspondence 
between the shift in the spectrum of optical modes 
and that of energy levels. Against this conjecture we 
have that there does not exist a priori a simple re- 
lation between the density of states for emission p3| . 



The approximate expression of <f>^^'^l for < 1 in Ea.l l96l l have 
been used for simphcity since the exponential factor in the inte- 
grand of Ea. l|95|l suppresses frequencies greater than c/(2^) and 
the singularities of the integrand locate around cjq <SC c/(2^) . 
The ^-dependent terms may be more relevant than the radiative 
ones if the host medium is highly dissipative within the frequency 
range of integration. 
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and that for the vacuum en- 
. In favor, we have aheady 



found that ^f^'' 

can be expressed in function of the elec- 
trical susceptibility only [Eg. ([5^ ]. Moreover, in the qua- 
sicrystalline approximation J-^^ is a function of Xj_j|((?) 
[Eqs. ([711175]) ] ■ Further, the qc approximation turns into 
exact in the long- wavelength of the effective medium the- 
ory of an MG dielectric. For this reason we will con- 
centrate on the MG model. Nonetheless, we must bear 
in mind that the MG model neglects recurrent scatter- 
ing in xmGj which might be relevant, even in the long- 
waveleiigth limit, in dense media for which S,'^ ^ p^^ - 
cf. [H, [sl] and Appendix [BJ 



A. The Lorentz-Lorenz shift 

The problem of the binding energy of an effective 
medium has been addressed by Bullough and Obada [8] in 
a molecular crystal and by ourselves in [23i] . In the latter 
reference it has been found that in the long-wavelength 
limit can be split into the energy of transverse bulk 
modes and that of LFFs and longitudinal bulk modes 
according to 



-pV _ -pV 
MG — Sch. 



Here, J-g^i^ is given in Eq. 



(115) 
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2^ 
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and 
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{paa)^eMG[ao] 



'mg "^ith the electro- 



(116) 

where xmg are functions of ao in the qc approximation. 
An identical expression has been reported by Bullough 
and Obada [1], who have interpreted AT^q as the elec- 
trostatic binding energy. However, the decomposition in 
Ea. (|115p . the identification of J^^^i^ with the radiative 
energy and the identification of AJ-^ 
static energy are all questionable. 

Regarding the radiative energy, it was found in Sec lVBI 
that not all the radiative en ergy is given by J^^^h ■ l'^ the 
first place, it was shown in [23] that >C_L(g) induces radi- 
ation on the surrounding dipoles around an emitter, act- 
ing as a mediator of non-radiative energy transfer. For 
this to be the case, C±{q) must appear coupled to G±{q) 
prior to integration in q. Clearly, the decomposition of 
Eqs. (|115p precludes this. On top of that, there is also en- 
ergy in -^il'aug carried by indirect radiation modes which 
are missing in -Fj^^ . It comes from the transverse modes 
within xwil) which are neglected in the long-wavelength 

limit [eg., the term — mentioned before Ea. (|103p 
belongs to indirect radiation]. 

Regarding the electrostatic energy, it is the EM vac- 
uum energy obtained in the electrostatic limit, c — > oo. 
Following [8], such a limit must be taken prior to the 
long-wavelength approximation, making this way radia- 
tive modes vanish. On the other hand, for q£^ the 



transverse susceptibility equals the longitudinal one be- 
cause only electrostatic modes enter x±{li — 0)- How- 
ever, as argued in the previous paragraph, C±{q^ = 0) 
has a physical meaning when coupled to G± acting as an 
inductor of radiation. Therefore, if radiation is precluded 
for c — > cxD, Q'^°° does not contribute to the electrostatic 
energy and neither do the transverse LFFs. As a result, 
transverse LFFs should not enter the vacuum energy in 
the electrostatic-long-wavelength approximation despite 
their presence in Ea. (|116p for q^ 0. Thus, we find that 
the energy of electrostatic modes in the long-wavelength 
limit is 
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which equals AJ^j^jQ but for the absence of the two trans- 
verse LFFs. 

Next, let us make explicit calculations in the MG 
model with bare polarizabilities. According to it, the 
electrical susceptibility is the sum of a geometrical se- 

(2) 

ries of ratio XmgI P'^^ — P^o/S, in which only the bare 
longitudinal propagator enters. Neglecting radiative cor- 
rections in Ea. (|10ip we get. 



XMG\Ola\ 



2/J_(£o 



P meo 



(118) 



The question we aim to address is whether the shift of 
the resonant frequency of xmg w.r.t. that of ao has any 
counterpart in the binding energy of the dipoles in the 
dielectric. That frequency shift is known as the Lorentz- 
Lorenz (LL) shift, and it is observable in optics experi- 
ments 55jp^. 



Aujll = ^^qU I - P 
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(119) 



This shift is clearly the result of the renormalization of 
Xmg in Ea. (lll8|) by the electrostatic modes within Xmg- 
Therefore, its energetic counterpart must find in the for- 
mula for J^Yff\stat- Since no radiative corrections enter 
the integrand of Eq. (jll7l) at all, it does not contain imag- 
inary terms and we can write it as a sum over modes. 

That is, upon using the regularization J (^2tt)'^ ~ ^ 
integral reduces to the sum of the poles minus the sum 
of the zeros of the factors within the logarithm. Using 
the formulas of an MG dielectric we obtain. 



^Yff\stat — P 




(120) 



Despite of the fact that there are experimental evidences for 
Aiii^^, it is also known that the MG model fails to provide ac- 
curate results close the resonance. The reason bein g th at it does 
not incorporate the dominant recurrent scattering [33(1 . 
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At leading order in p, J^Yfflstat — phAujLL holds. Be- 
cause, by assumption, no electrostatic-long-wavelengtli 
modes renormalize the single particle polarizability alone, 
we infer that J^Yfflstat is the binding energy of collec- 
tive d.o.f. due to long-range-electrostatic interactions. 
In a fluorescence experiment Aujll would be observed 
through the emission of several excited atoms of the same 
cluster decaying at a time. 

B. Discussion on the continuum approach 

Our microscopical approach together with the calcula- 
tions of the previous subsection allow us to answer two 
of the questions posed by BuUough in [i^l- The first 
one was whether the knowledge of the refractive index 
spectrum is sufficient to estimate the vacuum energy of 
long-wavelength fluctuations. Our results mean that nei- 
ther the radiative nor the electrostatic energy can be cor- 
rectly accounted for this way. The second question was 
concerned with the role of the Lorentz field in the contin- 
uum approach to the binding energy. We have proved, by 
strict isolation of the vacuum energy of long- wavelength 
longitudinal modes, that the Lorentz field gives rise to 
the LL shift and hence to J-'Yff\stat- 

Regarding the radiative energy, we conclude from the 
previous sections that the knowledge of the spectrum of 
the refractive index is insufficient for the quantification 
of the total radiative energy. In the first place, part of 
the radiation is not accounted for in the (bulk) spectrum 
of J-^g^i, ■ Second, the microscopical approach of SecfVl 
and Appendix shows that the transverse propagators 
within the LFFs terms amount to radiative modes which 
are disregarded in the long- wavelength limit of Eqs. (l73p - 

Regarding the electrostatic energy, we first observe 
that the relation J-Yff\stat — pfiAu^L must not be in- 
terpreted as a quantitative estimate of the binding en- 
ergy but as a result of consistency. That is, this equa- 
tion means that if only the longitudinal long-wavelength 
modes which enter the renormalization of xmg are con- 
sidered in the computation of the binding energy of clus- 
ters, the equivalence must hold. On the contrary, it 
must not be interpreted as a faithful estimation of the 
electrostatic binding energy. To see this it suffices to 
verify from Eg. ([TOP)) that pJ^o{p^)\hT'"'' does not con- 
tains terms of the order of Eq. (|120p . In particular, all 
the ^o-independent terms in that equation have a radia- 
tive origin and are of the order of Fq /wq less than that 
of TYff\stat- We conclude that, contrarily to what sug- 
gested in H^l, the continuum approximation is not 
even sufficient to estimate the orders of magnitude of 
electrostatic binding energies 



VIII. SUMMARY 

We have carried out a microscopical study of the EM 
vacuum energy of an isotropic and homogeneous molecu- 
lar dielectric made of a random distribution of two-level 
atomic dipoles. 

Before considering statistical averages, we have shown 
that for a specific configuration m, J-"^ can be expressed 
either as a function of the dipole fluctuations [Eg. ([5^ ]. as 
a function of the source EM field fiuctuations [Eg. ((53|) ]. 
or as a combination of both [Eq.((54|)]. 

When statistical averages are taken, the Lamb shift is 
a function of the electrical susceptibility only [Eg. ([55|) ]. 
Xj^ l|, and hence can be computed out of optical observa- 
tions. On the contrary, the total vacuum energy is not, 
and only a cluster expansion is possible [Eg. (p7)) ]. Only 
in the quasi-crystalline approximation it is possible to 
give a closed expression for the vacuum energy in terms 
Ofxr^ll [Eqs.dZlZSI]. 

Using a hard-sphere model we have discussed to which 
extent recurrent scattering terms contribute to the total 
vacuum energy. Except for the free-space Lamb energy, 
needed of a UV cut-off, no other divergences either in 
frequency or momentum space show up in the rest of 
the vacuum energy. In momentum space, LFFs kill the 
short distance divergences in the same manner they do 
in the spectrum of emission [23l . [53| . In frequency space, 
the UV divergences of retarded modes are exponentially 
suppressed by a natural cut-off at the wavelength of the 
order of the correlation length. On the contrary, both 
momentum and frequency divergences show up in the 
vacuum energy of an effective medium which is a func- 
tion of the effective bulk propagator only and contains no 
LFFs. For this reason we interpret that the LFFs play 
the role of spectral functions w.r.t. the spectrum of bulk 
modes in the sense introduced by Ford [56] . 

Related to the last issue, the Schwinger approach to 
the vacuum energy of an effective medium is in general 
insufficient to compute the Lamb shift out of variations of 
■^Sch ' which contradicts the result in [g^l -Secs lVBj and 
ED J^sch does not account for the total vacuum energy 
and contain artificial divergences because of not taking 
proper account of LFFs. Nevertheless, T^^f^ is sufficient 
to study Casimir forces between macroscopic dielectrics 

By evaluating the vacuum energy of the electrostatic- 
long-wavelength modes of a Maxwell-Garnett dielectric 
we have obtained that J-^fflstat — P^Awll- This rela- 
tion has been interpreted as a result of the consistency be- 
tween the spectrum of xmg and the energy of the modes 
involved in the construction of xmg- It has been con- 
cluded that the knowledge of the spectral form of the 



The sum over modes of Eq. llll6ll yields —phuJo[3 — 
2^-i/(nl + 2) - no^3/(ng +2)]/2 [I upon use of the regular- 



ization '' ..-i 



p and no = + XMci}^ = 0)- Equivalently, 
for a generic Lorentzian dielectric constant with field strength 
factor /, it yields pfnwo/2i 
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refractive index is insufBcient either to quantify the en- 
ergy of radiative modes or to estimate the electrostatic 
binding energy of a molecular dielectric. 
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Appendix A: Decomposition of cp^^''^^ and 
computation of <?f>o^'^' 



We write Eq. (|92|l as a spatial space integral in order to 
show the radiative and electrostatic nature of the prop- 
agators involved. To this aim, we make use of the fact 
that Cf^j^ (r) is totally transverse and G^^^j (r) is totally 
longitudinal in Fourier space, 



+ G^l.ir)-G^l.{r)hir-0 



dV 



(Al) 
(A2) 



+ Gi?LW-CWMr--0}, (A3) 



^i^) = -fcVaTrj J d^GZir) + G^.M 

■ Gil. W[l + Hr - 0]} = -k'pa ( d\ 



(0) 



^(0) 



X Tr 



{^G^^^t {r) ■ G[2n {r) 

+ GiZ.ir)-Gi'llXr)hir-0 
+ G[l(r)-GSLWM^-C)}- 



(A4) 

(A5) 
(A6) 



Note that it is the presence of h{r — ^) that makes the 
crossed terms non-vanishing. 

Next, particularizing to the hard-sphere model with 
h{r-0 = -e{r-0, 



~k 



pa[2 - 2iC + e'^ (-2 -|- 4iC + - iC^))] 



(1,0) I —'^ iC 



pa 



2n 
~k 
2^ 
-k 
2^ 



pa 



1 5 
2C ^ 6 
- 1 i 

1 1 i 



c<i, 



, C<i, 



(A7) 



(A8) 



where C = k^. The ^-independent terms correspond 
to the long-wavelength propagating modes. In Fourier 
space, they are given by the poles of the radiative prop- 
agators in Eas. (|AllA3IA6p . 

Next, we estimate the contribution of the two-point 
overdensity correlation function to ^^^'"^ in ^^fpy J'Q(^p2y 

Let h'""^ir - C) = (GS^^^^r - where G is a positive 
constant which accounts for the molecular coordination 
number. It suffices to integrate in angles to obtain. 



^(1,0) 
ova 



Tr 



d^rG^°\f){~k^pa)&"\f) 



2n 



■Cpae 



3 6 5 „ , 
^-^^-^ + 2* + C 



(A9) 



Appendix B: Computation of ^i^' and Xarec 
including recurrent scattering 



We proceed to sum up the infinite series of recur- 



^ofp)' -^o(p2)- "^^^ series is pictured diagrammatically 
in Fig|3^a). Replacing a with a the sum reads, 

0(^1) = _fc4^aTr{ j dVG(")(r) (Bl) 

• ^(fc2a)2™[G(0)(f)]2" •G(°'(f)[l + /i(r-0]}, 



where the scatterers are all taken bare to keep the order 
p in the series. An analytical expression for this sum can 
be given. We follow the computation of jSTl] in which the 
authors decompose G^^-* (r) in transverse and longitudinal 
components w.r.t. the position vector, r. 



G(°)(f) = P(r)[I-f(8)f] + Q(r)f®f, 



(B2) 



with 



P(r) 
Q{r) 



Airr 

_^tkr 

47rr 



[l + i/{kr) - l/{krf] 
[^2i/{kr) + 2/{krf]. 



In terms of P, Q, Eq. ljBll) reads. 



-k pa I dT 



2P2 



1 - (fc2aQ)2 



.1 - (Fq,P)2 
[1 + %-^]. 



(B3) 



The integration in r of Eq. (|Bip and its further integra- 
tion in u in Eg. (1911) require their expansion in powers of 
a. Hence, the series of Eq. ((98)) . 

In a similar fashion, the full series of recurrent scatter- 

(2) 

ing diagrams which amount to Xrec± \\ those on the 
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\ =h-h 



■■ 9=-k'pa 



• • • 



I I 



• • • 



(a) 



+ i 



FIG. 4. Diagrammatic representation of the series of recurrent 
scattering diagrams which amount to (a) (j>^^^ and (&) x'^'- In 
the latter series the first diagram corresponds to x'^'"' while 



,-,(2) 



the remaining ones belong to x 



right-hand side of the equahty of Fig|3|^6), but for the 



first one. Their sum yields, 

xi'L« = -fc'(P«)'G(o)(rO (B4) 



m—l 



Using the decomposition of G^"-' in terms of P, Q, it can 
be written as 



-(I - f f) 



1 - [k'^aPy 

[l + h{T-0]- (B5) 



1 - (fc2aQ)2 



At leading order in (a/^'^)^, the zero- mode which mod- 
ifies the original MG formula for the electrical suscepti- 
bility is 

X^^LW.M = 0) ^ i(pfi)2(5/4<3)2. (B6) 
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